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ABSTRACT 

Context. According to our current understanding, a subclass of the upper main sequence chemically peculiar stars, called mercury- 
manganese (HgMn), is non-magnetic. Nevertheless, chemical inhomogeneities were recently discovered on their surfaces. At the 
same time, no global magnetic fields stronger than 1-100 G are detected by modem studies. 

Aims. The goals of our study are to search for magnetic field in the HgMn binary system 66 Eri and to investigate chemical spots on 
the stellar surfaces of both components. 

Methods. Our analysis is based on high quality spectropolarimetric time- series observations obtained during 10 consecutive nights 
with the HARPSpol instrument at the ESO 3.6-m telescope. To increase the sensitivity of the magnetic field search we employed 
a least- squares deconvolution (LSD). We used spectral disentangling to measure radial velocities and study line profile variability. 
Chemical spot geometry was reconstructed using multi-line Doppler imaging. 

Results. We report a non-detection of magnetic field in 66 Eri, with error bars 10-24 G for the longitudinal field. Circular polarization 
profiles also do not indicate any signatures of complex surface magnetic fields. For a simple dipolar field configuration we estimated 
an upper limit of the polar field strength to be 60-70 G. For the HgMn component we found variability in spectral lines of Ti, Ba, 
Y, and Sr with the rotational period equal to the orbital one. The surface maps of these elements reconstructed with the Doppler 
imaging technique, show relative underabundance on the hemisphere facing the secondary component. The contrast of chemical 
inhomogeneities ranges from 0.4 for Ti to 0.8 for Ba. 
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1. Introduction 66 Eri (HR 1657, HD 32964, HIP 23794) is a well-known 

binary system with nearly i dentical components. U sing the ra- 

\^ The phenomenon of spot formation on the stellar surfaces is dial velocities determined bv lFrost & Struvel ([l924t) and his own 

. known for a long time. Detailed studies of stars exhibitmg sur- observations, dlSTl) found the orbital period of the sys- 

■ face features showed that the formation of temperature or chemi- he P ^ 5. 523 days, which was improved by the later 

Q ; cal spots is closely related to the presence of magnetic field in the studies (lYushchenko et al. 2001 ; Catanzaro & Let^ Ml). The 

^ stellar atmospheres (Donati & Landstreet 2009). Nevertheless, orbital analysis yie lds a mass ratio Mr /MaB very close to unity 

^ _ there are exceptional cases when spectropolarimetric studies of lYushchenko et al.ll200THc^anzaro & Letoll2QQ4l) . 

••. high sensitivity detect no magnetic field but spots are still present ^ ^. ^ ^ t- • . • i R7 1 

>► ^ n r T^u- • u J • u 11 1- Photometric studies of 66 En are controversial. lYoungl 

on the stellar surface. This is observed m chemically peculiar , . i .i r it t- i — 

n /r^r»\ ^ ^ /TT Ayr \^ (119761) has estimated the times of possible eclipses based on the 

(CP) stars or mercury-manganese (HgMn) type. ' , . ' , . r .1 . tt- 1 . . - 1 

rN r^, n . Txn./ , Orbital clcments 01 thc systcm. His photomctHC obscrvations at 

H The first discovery of sp2te in a HgMn star was reported ^^ ^^^^^ ^^^^ ^j^^^ variability. Young con- 

.9. for a And ^Adelman et alJ 2002|) which a so^ws an evo- ^j^^^^ ^j^^^ ^^^j^^l incUnation angle is less than 80°. Later, 

lution of spot structure ( |Kochukhov et al. | IZQQj. Currently ^^^^^^^ detected weak variability with a period dif- 

we know SIX other spotted H gMn starsLiiB_l If and ^^^^^^ ^^^^ ^^^j^^l one. This weak photometric variability 

HR 8723 ^Kochuldio^» AR Aur E^t al 200j ^^^^^^^ by Yush chenko et al. l who in turn sug- 

^Isom et al 2010), HD 1 1753, HD 53244 and HD 221507 variability Jth half of the orbital period. 

(IBriquet et al.i .20 IQ) . All these stars exhibit spots of mostly ^ ■ , . 

heavy chemical elements. According to theoretical models iBerghofer & Schmitt| (|1994|) reportedjhis star as an X-ray 

adopted for CP stars, physical processes responsible for the spot ^he later study | Hubrig et al. | (|200ll) detected the third, 

formation are associated with magnetic field. In order to enlarge low-mass companion, which is likely to be the actual source of 

the sample of spotted HgMn stars studied in detail, we have per- ^-'^y radiation coming from the system. This was directly 
formed a comprehensive analysis of another HgMn star. 



^ Throughout our study we denote the more massive component 

* Based on observations collected at the European Southern (HgMn star, the primary) with letter 'A". The less massive secondary 
Observatory, Chile (ESO program 084.D-0338). component is denoted by "B". 
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confirm ed with the Chandra X-ray imaging by IStelzer et"an 
(l2QQ6h . 

The early chemical analysis of 66 Eri by lYoungl ([1976|) 
showed overabundances of Hg, Cr and Y in the primary 
component, while the s econdary was reported to be normal. 
IWoolf & LambertI (1 19991) determined Hg abundance in a large 
number of HgMn stars, including 66 Eri. They confirmed the 
overabundance of mercury in the primary com ponent. A detailed 
abundance analysis of 66 Eri was presented by lYushchenko et al] 
(1999), who studied chemical composition of both components 
of the binary system. Combining spectroscopy and photometry, 
they determined Teff = 1 1 100 ± 100 K and T^^ = 10900 ± 100 K 
for the component B and A, respectively. The similarity of the 
physical parameters of the components also follows from their 
luminosity ratio Lb /La = 0.95 ± 0.05. The rotatio nal velocity 
of bo th stars is equal to Ve sin/= 17 kms"^ ( Yushch enko et al.l 
119991) . 

66 Eri is very similar to an other double-line binar y with an 
HgMn component - AR Aur dKhokhlova et aljri99^ . Despite 
the fact that no magnetic field was found in AR Aur, it shows 
chemical inhomogeneities (Folsom et al. 2010). The similarity 
in the two HgMn SB 2 stars inspired us to analyse 66 Eri with 
new high-precision observational data. 

The paper is structured as follows. In Sect. |2l we describe 
observational material and describe the data reduction process. 
Section [3] presents the least squares deconvolution analysis. In 
Sect. |4] we discuss the measurements of magnetic field in 66 Eri. 
Section [5] describes the spectral disentangling. Determination of 
the system parameters is given in Sect. |6] Sect. [7] presents an 
analysis of the spectral variability. Interpretation of the line vari- 
ability in terms of surface chemical maps is presented in Sect.[8l 
The paper ends with a summary of the main results and discus- 
sion in Sect. [9j 

2. Observations and data reduction 

Using the newly built polarimeter HARPSd oI (ISnik et al.ll2008L 
120101) attached to the HARPS spectrometer (iMavor et al.ii2003l) 
at the ESO 3.6-m telescope we have obtained high-quality obser- 
vations at ten phases, covering the full orbital period of 66 Eri. 
The resolution of spectra is ^= 115 000 and a typical S /N is 
200-300 at ^ ^ 5200 A. 

The HARPS detector is a mosaic of two 2K x 4K CCDs, 
allowing to record 45 and 26 polarimetric echelle orders on the 
blue and red CCD, respectively. The calibration set included 20 
bias exposures, 20 flat fields and two ThAr frames for each night. 
The flat fields and comparison spectra were recorded using the 
circular polarization mode. 

All science frames of 66 Eri were obtained in circular polar- 
ization, covering the wavelength range of 3780-6913 A with a 
small gap at 5259-5337 A. Each observation of the star was split 
into four sub-exposures corresponding to the 45°, 135°, 225° 
or 315° position of the quarter- wave plate relative to the beam- 
splitter. The length of individual sub-exposures varied between 
200 and 300 s, depending on the weather conditions. 

The ESO reduction pipeline was not available at that time , 
so we employed REDUCE package (Piskunov & Valentill2002h . 
written in IDL, to perform a number of standard steps to re- 
duce and calibrate our cross-dispersed echelle spectra. We av- 
erage and subtract bias images from the master flat and science 
frames. To locate spectral orders, REDUCE uses a cluster analy- 
sis method in conjunction with the master flat image. For the cor- 
rection of pixel-to-pixel sensitivity variations in the science im- 



ages the code normalizes the master flat field. Finally, REDUCE 
removes the scattered light and performs the optimal extraction 
of the stellar spectra as described by Piskunov & Valenti ( 20^. 

HARPS spectrometer is well-known for its stability, which 
is of the order of 1 m s"^ . Taking this into account, it is sufl&cient 
for our study to use only one ThAr spectrum obtained during 
each night to do the wavelength calibration. Employing ~ 700- 
900 ThAr lines we constructed a 2-D wavelength solution with 
an internal wavelength calibration accuracy of 18-21 m s"^ . The 
final step prior to the Stokes parameters derivation is continuum 
normalization. The continuum level is set by dividing each spec- 
trum by a smooth, slowly varying function, obtained by fitting 
the upper envelope of the blaze shape corrected, spliced spec- 
trum. 

For the derivatio n of Stokes V parame ter we use the ratio 
method described by iBagnulo et al.l (l2009h . This method mini- 
mizes spurious polarization eflTects by an appropriate combina- 
tion of the two physical beams recorded for four diflTerent re- 
tarder quarter- wave plate positions. Along with the polarization 
we also derive a diagnostic null spectrum. It is obtained by com- 
bining the polarization signal in individual sub-exposures de- 
structively, thus showing the residual instrumental polarization 
remaining after application of the ratio method. We use the null 
spectrum in the same analysis steps as the Stokes V spectrum, 
thus providing a realistic estimate of possible errors. 

Combination of eight spectra in four individual sub- 
exposures provided a convenient possibility to detect and remove 
remaining cosmic hits that otherwise seriously distort the final 
Stokes spectra. AflTected pixels are identified by their large devi- 
ation from the median and are substituted by the latter. Since our 
exposure times were relatively short, only 1-2 pixels required 
correction in each echelle order. 

3. Least-squares deconvolution 

Previous attempts to find magnetic fields in HgMn stars 
set an upper limit on the longitudinal magnetic field below 
100 G ('Shorlin et al."2002:'Wade et al.'2006:' AlI^iere et al.l201Ql : 
[Folsom etal. 2010; Makaganiuk et al. 2011). The fields of this 
strength are undetectable in individual spectral lines. In this 
case the sensitivity of magnetic field search can be increased 
by adding information from many sp ectral lines. Such a multi- 
line method has been implemented by lDonati et al.l (1199 7") in the 
so-called least- squares deconvolution (LSD). This technique as- 
sumes identical shape of spectral lines and represents observa- 
tions as a superposition of profiles scaled by appropriate weight. 
With this method it becomes possible to derive a high-precision 
average profile, increasing the signal-to-noise ratio by a factor 
of up to 1000. The latest studies of weak magnetic fields proved 
the LSD t echnique to be very eflTective in detecting fields weaker 
than 1 G (lAuriere et al.ll2009l; lu snieres et al. 1 I2QQ9 ). 

To construct a line mask required by LS D, we extracted 
the line list from VALD (Piskunov et al. 1995) using an atmo- 
spheric model with the eflTective temperature of 10900 K and 
abundance of chemical elements according to the results of 
lYushchenko etaP ([l999 ) for the more massive star. We note that 
LSD analysis is very weakly sensitive to the choice of Tq^ and 
abundances adopted for the mask. 

A line mask is a set of the laboratory line wavelengths and 
corresponding weights. In the case of the Stokes / spectra, the 
weight is a residual depth of individual line, given by VALD. 
For the Stokes V spectra the weight is the product of central 
depth, the line wavelength normalized by A = 4800 A, and an 
eflTective Lande factor of a given spectral line. Spectral lines with 
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Fig. 1. Variation of the LSD Stokes / profiles (left panel) and 
LSD Stokes V (right panel) with the orbital/rotational phase. For 
displaying purposes we shifted the Stokes / and V profiles verti- 
cally and multiplied the latter by a factor of 15. Phases are given 
on the left, above each profile. 



the central depth lower than a given cutoflT factor are excluded 
from the line mask. By setting the cutoflT criterion to 0.1, we 
produced the line mask contai ning 717 spectral lines. 

We use the LSD code by iKochukhov et alJ 12010) to com- 
pute LSD Stokes /, V, and null profiles. Profiles were computed 
in the velocity range from -170 to 200 kms"^ with a step of 
0.8 kms"^ This value corresponds to the average pixel scale of 
HARPSpol spectra. The errors of LSD profiles were calculated 
by propagating the uncertainties of individual pixels obtained at 
the reduction stage. 

In general there is a mismatch between a model LSD spec- 
trum and real observations due to coarse assumptions of the LSD 
technique. To account for these systematic errors we scale uncer- 
tainties of t he calculated LSD profiles by the reduced as de- 
scribed by Wade et al.l dlOOO"). The re- scaled LSD profile uncer- 
tainties are used for the estimation of the longitudinal magnetic 
field measurements. 

In Fig. [T] we present the LSD Stokes / and V obtained for ten 
orbital phases of 66 Eri. 



4. Magnetic field measurements 

The longitudinal magnetic field was determined from the firs t 
moment of the LSD Stokes V profile (IKochukhov et al.ll2Ql(ih . 
Due to the contribution of blends the continuum level of raw 
LSD Stokes / profiles is slightly oflTset from unity. We use a con- 
stant factor to correct the continuum level. This factor is applied 
to the LSD V and null profiles. Such a correction changes the 
value of the longitudinal magnetic field and its error by 2-3 G. 
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Fig. 2. Measurements of the longitudinal magnetic field with 
the corresponding error bars. Circles show the (5z) values for 
the primary component (HgMn star) and triangles represent the 
measurements for the secondary. 



For nine LSD V profiles we could measure the longitudi- 
nal field for each component separately. At the phase 0.827 (see 
Fig.[T]) we have measured both components as a single star. The 
results of the longitudinal magnetic field measurements are pre- 
sented in Table [T] Its columns contain the following informa- 
tion: the heliocentric Julian date of observation, the signal-to- 
noise ratio of the original spectra measured around 5200 A, the 
signal-to-noise ratio of the LSD V profiles, an orbital phase, de- 
termined from the binary solution discussed below. The next 
three columns represent measurements for the HgMn star (com- 
ponent A) as follows: the mean longitudinal magnetic field in- 
ferred from the LSD Stokes V, the longitudinal field inferred 
from the null profile, and a False Alarm Probability (FAP, see 
below). The last three columns give the same information for 
the component B. 

The errors of the magnetic field measurements were obtained 
by the standard error propagation, using the uncertainties of the 
LSD V profiles. The mean error of longitudinal field measure- 
ments is ~14 G for the HgMn component and ~20 G for the sec- 
ondary. For the HgMn star the longitudinal field value reaches 
the confidence level of 2cr only at the orbital phase 0.193. For 
the secondary star the measurements of magnetic field do not 
exceed -l.Scr level for most of the orbital phases. Thus, none of 
the individual (5z) measurements indicate the presence of mag- 
netic field in 66 Eri. 

Our magnetic field measurements for both components with 
their respective error bars are presented in Fig. 0as a function 
of orbital phase. The actual values of the longitudinal field mea- 
surements do not exceed the level of 3cr, so the periodic-like 
fluctuation of magnetic field measurements in both components 
are fully attributed to errors. 

Assuming a dipolar configuration of the magnetic field we 
made an assessment of the upper limit of its strength, which 
would be consistent with the entire set of (5z) measurements. 
Using a relation be tween the longitud inal field and parameters of 
an oblique dipole (iLeroy et al.ll 19941) . we estimated the strength 
of dipolar magnetic field to be 60-70 G with an uncertainty of 
50-70 G for both stars. However, based on the reduced statis- 
tics, the dipolar fit has no advantage over the null hypothesis. 
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Table 1. Magnetic field analysis of 66 Eri. 















A 










B 






HJD-24 X 10^ 


S/N 


S/N(hSD) 


Phase 


(5z>(V), 


G 


(5.>(null), 


G 


FAPxlO 


(5z>(V), 


G 


(5z>(null), 


G 


FAPxlO 


55202.725056 


245 


2638 


0.193 


28 + 


14 


7 + 


14 


4.92 


-26 + 


18 


39 + 


18 


5.49 


55203.744467 


299 


2674 


0.378 


30 + 


20 


-22 + 


20 


9.29 


16 + 


17 


8 + 


17 


7.73 


55204.737221 


320 


3304 


0.557 


-20 + 


17 


32 + 


17 


0.61 


2 + 


11 


-12 + 


11 


9.63 


55205.755838 


185 


2409 


0.741 


11 + 


23 


61 + 


23 


3.03 


-7 + 


16 


9 + 


16 


2.60 


55206.736354 


304 


2671 


0.920 


+ 


16 


15 + 


16 


2.59 


37 + 


23 


27 + 


23 


3.21 


55207.735135 


315 


2929 


0.100 


-2 + 


13 


-25 + 


13 


0.80 


32 + 


22 


-16 + 


22 


3.24 


55209.741479 


261 


2519 


0.463 


-10 + 


24 


3 + 


24 


9.75 


4 + 


15 


-13 + 


15 


1.76 


55210.746427 


216 


2225 


0.645 


-43 + 


26 


+ 


26 


2.11 


-19 + 


17 


5 + 


17 


9.74 


55211.750571 


241 


2965 


0.827 


8 + 


12 


+ 


12 


3.62 


8 + 


12 


+ 


12 


3.62 


55212.727307 


380 


3901 


0.004 


-1 + 


14 


40 + 


14 


8.28 


1 + 


11 


-14 + 


11 


6.58 



At the same time, much more complex magnetic fields can 
exist in the atmospheres of both stars. To investigate possi- 
ble presence of such complex magnetic structures, which might 
yield negligible (5z), we employed the EAP analysis of LSD 
profiles. EAP is a statistical estimate, which assesses if the 
observed structure in the Stokes V line profile is produced by a 
random noise. Following the definition by Don ati et aP (Il997h . 
these numbers should be interpreted as follows. No detection if 
EAP > 10"^ marginal detection if 10"^ <EAP< 10"^ and def- 
inite detection if EAP< 10"^. The smallest EAP obtained for 
HgMn star is 0.08. Thus, none of the measurements indicate the 
presence of a magnetic signal. A visual examination of each LSD 
Stokes y profile also do not suggest a magnetic field detection. 

To make sure that the LSD circular polarization profiles were 
not aflfected by a spurious polarization, we analysed LSD null 
spectra at each phase. The (^z) and EAP measurements were 
done in the same way as for LSD Stokes V profiles. The {B^) 
errors for the null spectra are fully consistent with those inferred 
from the circular polarization profiles. The absence of signal in 
the null spectrum means that our scientific data is not visibly 
affected by significant spurious polarization. 

5. Spectral disentangling 

Due to the orbital motion, the relative shift of the spectral lines 
of the components in the composite spectrum varies from zero 
to -200 kms"^ which can be seen from the LSD / profiles in 
Fig. [U In addition to this, the primary component shows an in- 
trinsic line profile variability. Since one of the main aims of this 
study is to investigate line profile variability in 66 Eri, it is im- 
portant to separate the eflTect of the spectral line blending due to 
the orbital motion from an intrinsic variability. 

To separate diflTerent variability effects we employ the proce- 
dure of spectral disentangling described by Eolsom et al. (2010). 
It combines information from all orbital phases and yields ra- 
dial velocities (RV) for both components, their average separated 
spectra and standard deviation in the reference frame of each 
component. Fig. [3] shows an example of spectral disentangling 
in the region 4355-4380 A containing several variable lines. 

The spectral disentangling procedure provided a new set of 
high-precision RV measurements for both components of 66 Eri. 
These data, presented in Table |2l have typical error bars of 
0.1-0.2 k ms~^ We coni b ined our RVs with the ve l ocities pub- 
lished bvlYound (Il976h lYushchenko etaP (1 19991 l2QQlh . and 
ICatanzaro & Letol (|2QQ4|) . deriving improved orbital parameters 
of the system. The new orbital elements are listed in Table [3l 
while Fig. 131 compares the respective predicted RV curves with 
the actual measurements. The standard deviation of the orbital fit 
is 2.5 kms"^ for the primary and 2.8 kms"^ for the secondary. 



Table 2. Radial velocities of the primary and secondary compo- 
nents of 66 Eri. 



A B 



HJD-24 X 10^ 


Phase 


yr,kms-^ 


Vr, kms"^ 


55202.725056 


0.193 


91.26 


+ 


0.19 


-27.92 


+ 


0.18 


55203.744467 


0.378 


-14.91 


+ 


0.20 


80.34 


+ 


0.14 


55204.737221 


0.557 


-60.30 


+ 


0.14 


127.11 


+ 


0.17 


55205.755838 


0.741 


-12.25 


+ 


0.13 


77.98 


+ 


0.16 


55206.736354 


0.920 


100.36 


+ 


0.11 


-37.29 


+ 


0.11 


55207.735135 


0.100 


132.44 


+ 


0.15 


-70.55 


+ 


0.13 


55209.741479 


0.463 


-46.91 


+ 


0.21 


113.26 


+ 


0.25 


55210.746427 


0.645 


-49.79 


+ 


0.09 


116.64 


+ 


0.16 


55211.750571 


0.827 


40.04 


+ 


0.09 


24.58 


+ 


0.21 


55212.727307 


0.004 


-74.02 


+ 


0.08 


135.94 


+ 


0.08 



Table 3. Results of the spectroscopic orbital analysis of 66 Eri. 



Parameter Value 



Porb (d) 


5.5226013+0.0000020 


7^0 (d) 


41356.499+0.017 


Kp^ (kms-^ 


102.83 +0.20 




100.35 +0.19 


y (km 


32.28 +0.10 


e 


0.0844 +0.0013 




160.9+1.1 




0.976+0.003 


MB/Mosin^"orb 


2.345+0.008 


Ma/Mq sin^ /orb 


2.403+0.008 


Ub/Rq sin /orb 


11.185+0.022 


aA/RQ sin /orb 


10.916+0.020 



Comp ared to the previous orbital analysis by ICatanzaro & Letol 
(I20Q4I) . we significantly improved the accuracy of the orbital pe- 
riod and of the velocity amplitudes, obtaining a mass ratio with 
the precision better than 1%. 

6. Fundamental parameters of the components 

The very similar properties of the 66 Eri components make it 
challenging to determine their parameters individu ally. The pre- 
vious detailed study of the system by Yushchenko et al.l (Il999|) 
suggested that the slightly more massive component, the HgMn 
star, is also more luminous, but cooler by 200 ± 140 K than the 
other star. The sign of this spectroscopically determined tem- 
perature diflTerence is impossible to reconcile with the main se- 
quence status of co-evolving binary components. 

We have redetermined the fun damental paramete rs of 66 Eri 
using Padova evolutionary tracks (iGirardi et al.l2000l) and taking 
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Fig. 3. Spectral disentangling procedure applied to 66 Eri. Dots show the observations and the solid line is the model spectrum for 
corresponding phase. Black solid lines at the bottom represent disentangled spectra for the primary and secondary components. 
The lines above the continuum level represent the standard deviations in the reference frames of components A and B. Peaks at 
the location of spectral lines in the standard deviation indicate the presence of intrinsic variability in a given line (Yii 4358.728 A, 
Ti II 4367.652 A, and Y ii 4374.935 A in this spectral interval). 



advantage of the improved orbital parameters determined in our 
study and an accurate Hipparcos parall ax, 7r= 10. 56 ± 0.34 mas, 
available for the system ( van Leeuwen & Eantino 2005). The 
primary component mass and the age of the system were found 
by interpolating in the set of evolutionary tracks for Z = 0.018 in 
order to match the total Hipparcos luminosity, L/Lq = 98.2+6.3, 
and the mean effect ive ter nperature, (Jeff) = 11000 K, deter- 
mined by Yushchen ko et al.l (11999) from the photometric prop- 
erties and spectral energy distribution of 66 Eri. The outcome of 
this analysis is illustrated in the upper panel of Fig. [5] Tabled 
summarizes the resulting fundamental parameters of the com- 
ponents, with error bars reflecting the uncertainty of the total 
luminosity. Efl'ective temperatures of the components are given 
without error bars because they are mai nly sensitive to (r eff), fo r 
which no uncertainty was provided by Yushchenko et al.l (1 19991) . 

Our analysis of the H-R diagram position of the components 
of 66 Eri suggests that the primary is more massive, hotter, and 
more luminous (La/ Lb = 1.093), as one would expect for two 
co-evolving main sequence stars. Our temperature difl'erence, 
Te if jA) - Teif(B) = 163 K, is opposite to what was suggested 
by lYushchenko et al.l (Il999l) . Nevertheless, the results of their 



abundance analysis of 66 Eri components are unlikely to change 
substantially due to ~ 200 K change of Tqq. 



The evolutionary masses and radii allow us to determine 
the orbital and rotational inclination angles, testing synchroniza- 
tion of the system. We found the projected rotational veloci- 
ties, Ve sin /rot = 17.1 ± 0.2 and 16.9 ± 0.2 kms"^ for the com- 
ponent A and B, respectively, from fitting the profiles of ~20 
unblended lines in the 4400-4600 A wavelength region. These 
results are in excellent agreement with Ve sin /j-ot — 

17 kms"^ 

given by Yushchenko et al. (1999) for both stars. Since the radii 
are known, the oblique rotator relation yields /rot ~ 74° assum- 
ing Prot = ^orb = 5^5226 (scc lowcr panel in Fig. O. On the 
other hand, from the orbital solution and evolutionary masses 
/orb = 76°. This convergence to nearly identical inclination an- 
gles derived using independent methods is a strong indication 
that the rotation of the 66 Eri components and their orbital mo- 
tion are coaxial and synchronized. 
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Fig. 5. Luminosity-effective temperature diagram constructed 
for the components of 66 E ri with the evolutionary calcula- 
tions by iGirardi et al.l (l2000b . The thin solid curves represent 
the original evolutionary tracks, while thick curves show the 
tracks interpolated for the final component masses. The dotted 
line shows the isochrone for log t-l Al yr. The vertical dashed 
line corresponds to mean effective temperature of the compo- 
nents, Teff = 11000^. 

7. Line profile variability 

The line profile variability search was performed based on the 
disentangled spectra of both components. We carefully exam- 
ined each spectral region, excluding hydrogen lines and the re- 



Table 4. Fundamental parameters of the primary and secondary 
components of 66 Eri. 



Parameter 


A 


B 


log t (yr) 


7.47 + 0.23 


MIMq 


2.629 + 0.032 


2.566 + 0.032 


R/Rq 


1.948 + 0.063 


1.919 + 0.061 


L/Lq 


51.3 + 3.3 


46.9 + 3.0 


TefF (K) 


11077 


10914 


logg 


4.28 + 0.03 


4.28 + 0.03 


Ve sin/ (km s~^) 


17.1 +0.2 


16.9 + 0.2 



gions populated by the telluric lines. To diagnose the variability 
in spectral lines we computed standard deviation in the refer- 
ence frame of each component. It is shown in Fig. [3] by a line 
above the average spectra of components A and B. An excess of 
standard deviation at the positions of spectral lines indicate their 
variability. 

The variability search in spectral lines was performed for 
both components. We confirm that the secondary component 
does not exhibit any variability, while for the primary star we 
found variability in the spectral lines of Ti ii, Sr ii, Y ii, and Ba ii. 
This variability is generally weaker than found in other spotted 
HgMn stars. We suspect that a weak variability is also present in 
some Cr lines. 

Earlier st udies of HgMn stars (l Adelman et al.l I2QQ2I; 
iKochukhov e t al. 2005) focused on the variability of Hg ii line 
at /I 3984 A. Our analysis of this spectral line shows that, de- 
spite being rather strong, it shows extremely weak variability in 
66 Eri A. 

The full list of variable lines is given in Table[5] The first col- 
umn of this table gives a line identification, the second lists the 
corresponding laboratory wavelength. The last two columns pro- 
vide the same information for occasional blending lines. Lines 
with marginal variability are marked with italics. 

In Fig. Owe show typical variable spectral lines of the four 
chemical elements with definite line profile changes. Spectra 
are plotted together with an average profile (dotted line). The 
strongest variability is found for Sr and Ba lines; the weakest is 
seen for Ti. Examination of the line profile variability showed 
that all spectral lines listed in the Table [5] show the largest pro- 
file distortion relative to the mean spectrum at phases 0.193 and 
0.378. The variability of Ban is also clearly seen for the phase 
0.463. In addition, the lines of Yii and Srii also marginally de- 
viate from the corresponding average profiles at phases 0.827, 
0.920 and 0.741, 0.827 respectively. 

The observed spectral variability occurs with the period of 
the orbital motion. For all four chemical elements with variable 
lines we found a distinct change of the radial velocity and, es- 
pecially, equivalent width with the orbital period. Since our ob- 
servations cover only two rotation cycles, we cannot determine 
variability period with a high precision. However, we observed 
that the integral properties of variable lines and their shapes re- 
peat closely after one orbital period. This strongly suggests that 
the orbital motion and the stellar rotation are synchronized. 

8. Doppler imaging 

To study surface distribution of chemical elements in 66 Eri 
A, we applied the Doppler imaging (DI) method. We used 
INVERSIO code (Piskunov & Kochukhov 2002) to determine 
abundance maps, ignoring magnetic field. This analysis was ap- 
plied to several non-blended spectral lines of Ti and Y that ex- 
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Fig. 6. Illustration of the line profile variability in the primary component (HgMn star) of 66 Eri. On the panels from left to right we 
show spectral lines as follows: Sr ii 4215. 159 A, Ti ii 4779.985 A, Ba ii 4934.076 A, and Y ii 5200.410 A. The dashed line represents 
an average profile. All profiles are shifted upwards for displaying purposes. 



hibit a significant variability and to all available spectral lines of 
Sr and Ba. In total, we used two spectral lines of Srii, two of 
Ban, six lines of Y ii, and nine spectral lines of Ti ii. These lines 
are marked with bold font in Table |5] Multi-line DI is applied 
here to study spots in an HgMn star for the first time. This al- 
lowed us to reconstruct the surface maps with a higher precision 
than what was usually achieved using a single spectral line. 

Doppler inversions were performed using the inclination an- 
gle of 75°, as determined earlier (Sect. O. The initial abun- 
dances of Ti, Sr, Y, and Ba were taken from lYushchenko et al.l 
(Il999l) . Spectra were calculated employing the revised model 
atmosphere of 66 Eri A with parameters Teff = 11100 K and 
logg = 4.25. The rotational velocity was adjusted to give a bet- 
ter fit to the observed line profiles, yielding the final value of 
Ve sin / = 17.5 ± 0.2 km s"^ . The reconstructed surface maps are 
shown in Fig.[7l Online Figs. [81-[TQ| illustrate the comparison of 
observed and model line profiles. 

First of all, we can see from these maps that there is a gen- 
eral asymmetry of chemical composition for two hemispheres of 
the star. For all four studied chemical elements we find a larger 
abundance at rotational phases from 0.5 to 0.1 and a smaller 
abundance for phases 0.1-0.5. Secondly, there is a diflTerent spot 
structure for the hemisphere with element overabundances. 

The spots of Ba and other chemical elements except Ti ex- 
tend from -30° to 60° in latitude, showing the highest concentra- 
tions at the phase 0.65. For Y and Sr we also find the secondary 
spot at the phases 0.95-1.05. The spots of Sr and Ti are more 
extended azimuthally, while the spots of Ba and Y are stretched 
meridionally. Sr is also the only element of the four that shows a 
relative overabundance at the pole. Finally, Ti shows the spot po- 
sitions similar to Sr, but with a discontinuity at the stellar equator 
even for the region of relative underabundance. 

For all chemical elements presented in Fig. [7] the abundance 
diflference between zones of the highest and lowest concentra- 



tions does not exceed 1 dex. For Ba the range is 0.81 dex, for Y 
it is 0.71 dex. For Ti and Sr the range is 0.36 dex and 0.59 dex, 
respectively. This is in a good agreement with the relative am- 
plitude of variability in individual line profiles of these chemical 
elements (Fig.O. 

To exclude the possibility that spots are found by DI code at 
roughly one hemisphere of the star due to a spurious systematic 
eflfect, we applied the same DI analysis to ten spectral lines of 
Fe, which show no variability. The inferred abundance map of 
Fe is not similar to the maps of other chemical elements in its 
morphology. In particular, it does not show an abundance dif- 
ference between two hemispheres. The diflTerence between the 
regions of the highest and lowest abundances does not exceed 
^ 0.15 dex, which we take to be an upper limit of the possi- 
ble inhomogeneity for chemical elements without obvious line 
profile variability. 



9. Summary and discussion 

Earlier studies of several HgMn stars found a variability in some 
spectral lines, attributing it to the presence of chemical spots. 
However, all previous attempts to find magnetic fields that can be 
responsib le for these chemical inhoni ogeneities we re unsiiccess- 
ful (ShorHnetaP l2002t IWade et al.l .2006: Fols om et al.l 120101: 
lAuriere et al.ll2010l: iMakaganiuk et al.ll201 ih . 

To extend the group of spotted HgMn stars with known con- 
straints on magnetic field, we performed time-series spectropo- 
larimetric observations of 66 Eri. This is the first such analysis 
of this unique spectroscopic binary star. It allowed us to search 
for magnetic fields in both components and investigate their line 
profile variability. 

Taking the advantage of one of the best spectropolarimeters 
in the world, we recorded high S /N and high-resolution spectra 
of 66 Eri covering its full orbital period. Despite the high polari- 
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Fig. 7. Surface distributions of Ti, Sr, Y, and Ba abundances reconstructed for 66 Eri A. The star is shown at five equidistant rotation 
phases at the inclination angle / = 75° and vertically oriented axis of rotation. The darker spots represent a higher concentrations 
of chemical elements. The contours of equal abundance are plotted with a step of 0.1 dex for all chemical species. The thick line 
indicates the position of rotational equator, while the vertical bar shows rotational pole. The scale of abundance maps is shown on 
the right. Abundances are given in a logarithmic scale relative to the Sun. 



metric sensitivity of the instrument, we did not find any evidence 
of the longitudinal magnetic field in either component of 66 Eri. 
Our measurements of (5z) using LSD profiles have error bars 
10-24 G for both sta rs. This null result agrees with the outcome 
of our earlier work (iMakaganiuk et al. where the longi- 

tudinal magnetic field was measured in a sample of 47 HgMn 
stars. The analysis of these objects showed no evidence of weak 
magnetic field signatures in Stokes V spectra. 

Our magnetic field measurements provide an upper limit of 
60-70 G for a dipolar field in either star and LSD V profiles pro- 
vide no evidence for the presence of strong, complex magnetic 
field structures. This upper limit is significantly smalle r than 
the rn inimum dipolar field diagnosed in Ap stars (Auriere et al.l 
[2QQ7h . Thus, 66 Eri is clearly not an Ap star. 

As we deal with the binary system that shows two systems 
of spectral lines in its composite spectrum, it is necessary to dis- 
entangle the spectra to study each component separately. This 
was accomplished with the procedure of spectral disentangling. 
Along with the intrinsic spectra of the components it provided 
their radial velocities. The fit of our RV measurements together 
with RVs from previous studies (lYoungl 19761: lYushchenko et al.l 



ll999Ll200ll: ICatanzaro & LetQll2004 with theoretical model, en- 
abled us to improve the orbital elements of 66 Eri. 

Using constraints from the orbital solution, we found the fun- 
damental parameters for each star with the help of evolutionary 
tracks. The primary (HgMn) star turned to be more luminous and 
hotter, as expected. Yushchenko et alj (1 19991) considered HgMn 
star to be the secondary component in their abundance analy- 
sis. The eflTective temperature revision of 200 K suggested in our 
study is relatively unimportant for abundance determination. 

The evolutionary tracks provided masses of the components, 
which were utilized to compute the orbital inclination angle. The 
other method, based on the oblique rotator model, yielded nearly 
the same value for the rotational inclination angles. These results 
indicate that the angle between the rotation axes of two stars and 
the orbital plane is 90° and their rotation is synchronized with 
the orbital motion. 

For the identification of spectral lines, we computed a 
synthetic spectruni b ased on the abundances provided by 
Yushche nko et al.l ('1999) for both components. For the primary 
component the agreement between our synthetic spectrum and 
observations is rather good in most cases, in particular for the 
lines of Fe, Ti, Cr, Y, Ba, Si, Mn, and Hg. At the same time, we 
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Table 5. The list of variable lines in the primary component of 
66 Eri. Bold font marks spectral lines used for Doppler imaging 
and italics indicate lines with marginal variability. 
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cannot confirm th e high abundances of La , Zr, Yb, Hf, Pt, Au, 
and W inferred by lYushchenko et all (1 19991) . The strongest lines 
of these elements are very weak or absent in the observations 
while in synthetic spectrum they are too strong. Th e abundances 
of all these elements, derived by lYushchenko et aD (1 19991) based 
on a small number of very weak spectral lines, are not reliable. 



The abundances for the secondary component seem to be in a 
very good agreement with our observations. 

In conclusion of their abundance analysis lYushchenko et al.l 
(Il999 ) made a confusing suggestion that the more massive 
component of 66 Eri is not an HgMn star. This claim con- 
tradicts the standard definition of an HgMn star as a late-B, 
slowly rotating ob ject wit h enhanced Hg or Mn lines (e.g. 
|Preston lll974.T^ shchenko et al.1 (Il999h themselves found both 
Hg and Mn to be overabundant; a l arge H g overabundance 
was confirmed by IWoolf & LambertI (1 19991) . The absence of 
P overabundance is unusual, but is kn own for several othe r 
HgMn stars with Teff similar to 66 Eri A (lAdelman et al.ll200lb . 
while a significant Ga overabundance is mostly inferred from 
the optical spectra of HgMn stars much hotter than 66 Eri A 
fey abchikova et al.l IT996) . Thus, 66 Eri A is undoubtedly an 
HgM n star and was treated as such in many modem studies 
(Cowley et al. 2007; Catanzaro & Leto 2004"; iRottler & YoungI 
[2002:.Woolf& Lambert 1999: Hubrig 1998). 66 Eri is also clas- 
sified as "Hg" in the updated "Ca talogue of Ap, HgMn and Am 
stars" (iRenson & Manfroidll2QQ9l) . 

The analysis of spectral lines in the disentangled spectra of 
66 Eri A allowed us to discover the variability in some of them. 
This was not achieved in the previous studies of this star because 
of a low resolving power and poor S /N. On the other hand, same 
analysis revealed no variable spectral lines for 66 Eri B. 

Spectroscopic variability of 66 Eri A turned out to be diff'er- 
ent from other HgMn stars studied previously with time- series 
spectroscopy. All of them show not able variability in Hg ii spec- 
tral lines. For example, for a And iKochukhov et"aD (^007) re- 
ported ^ 30% equivalent width changes of the Hg ii 3984 A line. 
The variability in this line is barely detectable in 66 Eri A. Our 
spectra are also insufficient to unambiguously show variability in 
Cr lines, although we found marginal changes in four of them. 

More pronounced variability is found for the spectral lines 
of Y, Sr, Ba, and Ti. In Fig. [6] we illustrate the line profiles of 
these four chemical elements that possess the strongest variabil- 
ity among all spectral lines of each species. The change in pro- 
files with the rotational phase is still rather weak compared to 
other spotted HgMn stars, which makes this star the least vari- 
able. 

We found slightly different line profile variability pattern for 
four elements. This indicates diff'erence in the distribution of 
these chemical elements over the surface of 66 Eri A. We applied 
Doppler imaging technique to study spot geometry. The derived 
surface maps (Fig.jT]) indicate higher concentrations of chemical 
elements on one hemisphere of the star. The overall abundance 
contrast is about 0.5 dex for all studied chemical elements. The 
hemisphere with a lower abundance still has higher concentra- 
tions of these elements compared to the Sun. This means that a 
certain physical process leads to a systematic difference in abun- 
dances between the two hemispheres. This phenomenon cannot 
be attributed to the strong magnetic field due to its absence in 
the atmosphere of 66 Eri A. 

Low-contrast abundance inhomogeneities found on the sur- 
face of 66 Eri A may give rise to photometric variability. So 
far the evide nce of li ght variation in 66 Eri was controver- 
sial. Yushchen ko et al. l (1999) could not verify weak photomet- 
ric variability suggested by (Schneided ll987h . but suggested the 
presence of variability with half of the orbital period in the 
Hipparcos epoch photometry of 66 Eri. Our reanalysis of these 
data does not support the claim by Yushchenko et al. There 
is a peak in the amplitude spectrum at P = 2.754 ± 0.004 d, 
which is close to, but formally is not consistent with, Porb/2 = 
2.761 d. However, this periodicity is not statistically significant 
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(FAP = 0.4) and there are more than ten other peaks in the am- 
plitude spectrum with the same or higher magnitude. 

Since 66 Eri is a binary system, we cannot exclude the in- 
fluence of the secondary component on the atmosphere of the 
primary. To understand how the position of chemical spots in 
66 Eri A correlates with the relative orbital motion of the com- 
ponents, we combined the results of the orbital analysis and DI 
maps of yttrium. The resulting picture is shown in Fig. [TT] 

It turns out that the relative underabundance of Ti, Sr, Y, and 
B a in 66 Eri A is observed on the hemisphere which faces 66 Eri 
B. Due to a rather close distance between the two stars, the tidal 
eflTects can be responsible for the destruction of the upper layer 
containing the excess of studied elements. On the other hand, 
the radiation flux from 66 Eri B can also have an impact on the 
distribution of chemical elements in the atmosphere of the pri- 
mary component. However, it is not clear why is only one of the 
nearly identical stars afl'ected by these efl'ects. 

In addition to the puzzling spectroscopic variability in 
66 Eri A, the two stars have significantly diflTerent average chem- 
ical characteristics. Despite the two stars having nearly the same 
fundamental characteristics, the secondary component seems to 
be a normal star with nearly solar chemical composition and 
shows no surface inhomogeneities. Remarkably, the line profile 
variability in A component occurs only for those elements which 
show much stronger spectral lines in the primary. This is also 
cle arly seen from the abunda nces of Ti, Sr, Y, and Ba determined 
by lYushchenko et aTl (1 19991) : all these four elements are signifi- 
cantly more abundant in the atmosphere of the primary. Among 
the elements with the largest abundance diflference between the 
components only Hg does not clearly show a non-uniform dis- 
tribution. 

The identical age of both stars and the small separation be- 
tween them suggest that they were formed in the same proto- 
stellar gas cloud. This would imply the same overall chemi- 
cal composition for both stars. The observed chemical inhomo- 
geneities in one of the two components and a large diflTerence 
between their average abundance could be related. These eflTects 
cannot be attributed to the generic atomic diflTusion, which de- 
pends only on fundamental stellar parameters. It must be a result 
of some rapid process that changes the atomic diflTusion on a 
timescale much shorter than the age of a star. If this is true, the 
secondary component has either already undergone the stage of 
non-uniform surface distribution of chemical elements or may 
show surface structures in the future. 
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Fig. 11. Orbital motion and rotation of the components in the binary system 66 Eri. The star with yttrium spots represents 66 Eri A 
and the other filled with a solid color represents 66 Eri B. The radii and the separation between the two stars are shown to scale. 
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Fig. 8. Comparison of the observed (symbols) and model (solid lines) profiles for Sr and Ba lines in 66 Eri A. The bars in the lower 
left corner indicate the horizontal and vertical scales (0. 1 A and 5% of the continuum intensity respectively). 
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Fig. 10. Same as Fig.[8]but for Ti. 



